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Abstract 

We consider the static and spherically symmetric field equations of 
general relativity for charged perfect fluid spheres in the presence of a 
cosmological constant. Following work by Florides (1983) we find new 
exact solutions of the field equations, and discuss their mass radius ratios. 
These solutions, for instance, require the charged Nariai metric to be the 
vacuum part of the spacetime. We also find charged generalizations of the 
Einstein static universe and speculate that the smallness problem of the 
cosmological constant might become less problematic if charge is taken 
into account. 



1 Introduction 

Ever since Schwarzschild constructed the first exact interior solution of the Ein- 
stein field equations, static and spherically symmetric matter solutions have 
been of great interest. Additional degrees of freedom of the matter were in- 
cluded at later stages, like the cosmological constant, anisotropic pressure and 
charge. An entire new class of interior Schwarzschild like solutions was obtained 
by Florides [11], these particular solutions were distinct since they featured a 
vanishing radial stress TJT. This Schwarzschild like solution was later general- 
ized to include a non-zero charge distribution, see jT2][T31[TB] . In each case the 
new interior solution was matched to an exterior Reissner-Nordstrom solution, 
or the exterior Schwarzschild solution in the absence of charge. 

Static charged fluid spheres had been studied in detail by Efinger (1965), 
for instance, who was interested in understanding charged particles in general 
relativity, see also [21 [HI IH1 HH1 H31 [151 US] - Note that the spacetimes obtained 
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so far were not singularity free. Then Mehra (1979) managed to construct an 
interior solution and avoid a singularity. Florides (1983) also considered charged 
perfect fluid solutions, these solutions described the complete field, interior and 
exterior. The neutral Florides interior solution was then extended further by [20) 
to include the effects of a non zero cosmological constant A. This solution was 
matched to the exterior Schwarzschild-de Sitter (Kottler) solution. Since then 
not much attention has been devoted to analyze charged perfect solutions in the 
presence of a cosmological constant. 

In this paper we follow |llj and impose the restriction that the fluid com- 
ponent of the matter is isotropic. In this way we obtain new classes of interior 
charged solutions in the presence of a cosmological constant. A particular class 
of these solutions, namely the charged analogue of the interior Nariai solution 
is discussed in greater detail. We show how to match this solution to the exte- 
rior charged Nariai solution. Note that the exterior is not asymptotically flat 
and thus has a somewhat special place in the family of Reissner-Nordstrom 
de Sitter solutions. We are also able to construct charged generalizations of 
the Einstein static universe where the pressure is monotonically decreasing or 
increasing between the two regular centers. 



2 Field equations 

We consider the system of Einstein-Maxwell equations with cosmological term 
Gij + Kgij = 8nTij , Ty- = My + , 

dj ( V^gF^ ) = j—g r , d {i F jk] = , 

where we set G = c = 1. The matter energy- momentum tensor is denoted by 
My , and the electromagnetic energy-momentum tensor Eij is given by 

E - = — ( F- k F u - -a F mn F 

^ij — ^ \ ± i ± jk ^JtJ ± ran 

The 4-current density vector J % is defined by 

J 1 = Airau 1 . 

The proper charge density is denoted by a and u 1 is the 4-velocity satisfying 
U l Ui = —1. 

Following the notation of [13] . we consider the static and spherically metric 
in the form 

ds 2 = -e a ^dt 2 + e b ^dr 2 + r 2 dtt 2 , (2.1) 

where dSl 2 = d8 2 + sm 2 9 deft 2 is the line element of the unit sphere. 

Using this metric ansatz and the above defined 4-current density vector, the 
only non- vanishing component of the electromagnetic tensor Fij becomes 

p (a+b)/2 rr 

F rt = ~^^<1 , q = 4tt / ae h ' 2 r' 2 dr' . (2.2) 
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Therefore, the resulting components of Ef are 



E l = K = -El = ~Ef = -£r 4 . (2.3) 

As for the matter we choose an anisotropic perfect fluid 

M? = diag(-p,p,px,px) . (2.4) 

The resulting (tt) and (rr) components of the Einstein-Maxwell field equa- 
tions and the energy-momentum conservation equation take the following form 



^(r--- b )-A = 8 7 rp+^, (2.5) 
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e-° a'e-° 1 , q 

ft + + p) + % - p±) - r^i = • ( 2 - 7 ) 

Therefore, this system consists of three independent equations for six unknown 
functions {a, 6, p,p,p±,q}- Therefore the system is under-determined and addi- 
tional physical or mathematical assumptions need to be taken into account in 
order to close the set of equations. 

For example, in the neutral case q = 0, with an isotropic perfect fluid p = p± 
there are three equations and four unknowns. In that case the most physical 
approach is to consider a barotropic equation of state for the matter, p — p{p) 
which closes the set of equations. Another popular approach motivated by astro- 
physical observations is the prescription of a density profile p(r). However, this 
usually leads to a singular center and is therefore problematic from a rigorous 
point of view. 

In the charged case, one can impose various physical conditions to decrease 
the number of unknown functions. One can for instance assume that the total 
energy momentum tensor is isotropic 



which yields the condition 



Tr = T$ = Tt , (2.6 



Therefore, this assumption leads to a relation between the anisotropy (difference 
between radial and tangential pressure) and the electromagnetic field. Note that 
in the absence of charges this immediately leads to an isotropic perfect fluid. 
The second natural assumption is to consider an isotropic perfect fluid 

M r r = M e e - Aft , (2.10) 
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which simply gives 

P = P±. (2.H) 

Note that more conditions need to be imposed to close the system of equations. 

A useful linear combination is the sum of equations (|2.5p and (|2.6p which 
yields 

a' + b' = 8Trre b (p + p) . (2.12) 
The first field equation (|2.5p can be formally integrated and we obtain 



2m q 2 A 9 
1 - + 2.13 

r r A 3 

where we used (12.21) and where 



rn 



+ m q = / 4irpr' 2 dr' + / inr' aqe b/2 dr' . (2.14) 
Jo Jo 

The field equations can be combined into the useful Tolman-Oppenheimer- 
Volkoff equation, which is a differential equation for the pressure expressed in 
terms of the other matter quantities. Eliminating a' from (|2.7[) using (|2.6p 
and (j2~T3)) . we find 

P ^ r i_^ + 4_A r2 r yp p^) 4^4- i 2 - 15 j 

- 1 - r ~ 7 1 3 

which reduces to the well-known Tolman-Oppenheimer-Volkoff equation when 
q = 0, A = and p = p±. 

If the interior solution is of finite extent, we denote total quantities with 
capital letters, this means M, Q denote the total mass and total charge of the 
solution, respectively. We use r b to denote the boundary of the object, defined 
as the vanishing pressure surface. Therefore, if the exterior spacetime is the 
Reissner-Nordstrom de Sitter spacetime, its metric is given by 

r r z 3 / 

i _2M + Q2_A r2 \ 1 dr 2 + r 2 dn 2^ (216) 
r r z 3 / 

where r > r^. 

The main aim of this work is to find new exact solutions of the field equations. 
Therefore, in view of Eq. (|2.13l) one particular choice would be to consider the 
case when 

q 2 

8np H — j = 8nfj, = constant . (2-17) 
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This condition has been first considered in [3] where a charged dust was con- 
sidered. It has further been exploited by [13] , In this case the interior metric 
takes the simpler form 



e-» = l- w , (2.18) 
1 Snn A 

Combining equations (12.19)) and (|2.18l) with (|2.17p yields the metric function 
e~ b in the form 

3 Solutions of the field equations 
3.1 Florides' solution 

Let us consider an isotropic total energy- momentum tensor (Condition 1) which 
gives (12.91) . Inserting this into the conservation equation (|2.7[) gives 



p' + ~A P+ p) + ^-±% = o. (3.1) 

2 2ir r° An r 4 



Now we substitute a' using (|2.12[) and use the fact that 

q 2 \ 1 q 2 1 qq' 



dr V Snr 4 J 2tt r 5 Ait r 4 
we arrive at the equation 

P > \b'{ P + P ) + W(p + vf I = o . 

Finally, we introduce the quantity z = p + p which results in 

z' -\b'z + A-Kre b z 2 --^-(p+^—]=Q. (3.2) 
2 dr \ 87rr 4 / 

If we assume the term in the bracket to be a constant (Condition 2), sec 
the condition (|2.17l) . then this equation simplifies considerably which we will 
assume henceforth, which yields 

z' - h)' z + Anre b z 2 = 0. (3.3) 

Next, making use of (|2.18[) we have 



2r/R 2 , v 
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and therefore 



r/R 2 



1 - r 2 /R 2 



1 



l /R- 



■z 2 = 0. 



which can be integrated using separation of variables 



4vri? 2 



1 + 



where C is a constant of integration. 

The other metric function a follows from (|2.12l) 



o«/2 _ 



l + CJl 



1 + C* 



(3.5) 



(3.6) 



(3.7) 



where we fixed the other constant of integration by requiring e a (r = 0) = 1. 
Therefore we found the complete metric. Note that in a static spacetime we can 
always rescale the time coordinate arbitrarily. 



3.2 Towards new solutions 

In [13] , the constant of integration C was determined by the condition that the 
pressure vanishes at some radius p{rb) = 0. Since this condition requires the 
existence of a vanishing pressure hypersurface, we follow a different approach 
and fix the constant of integration using the central values of pressure and energy 
density. 

Let us define the central pressure and energy density to be p c = p(r = 0) 
and p c = p(r = 0), respectively. Then 

Pc+Pc = z c (r 

C = 

In the chosen coordinate system there is a coordinate singularity when r — > 
R. The spatial part of the metric is a 3-sphere, hence it is natural to introduce 
the third Euler angle of the sphere by 

r 2 

r = Rsinx, 1 — = cos 2 x ■ (3.10) 

In these new coordinates the interior metric takes the form 

ds 2 = - ( 1 + ° C ° S X ) 2 dt 2 + R 2 d X 2 + R 2 sinVO 2 , (3.11) 

which covers the entire 3-sphere. Note that the matter part of the solution 
does not necessarily have to 'fill' this three sphere. In fact, solutions analogous 



) = Wi + cr (3 - 8) 

*iirR z {p c +p c ) 
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to the Schwarzschild interior solution have their vanishing pressure surface lo- 
cated before the equator of this three sphere and thus this new coordinate is 
not important in understanding them. However, in the presence of a cosmolog- 
ical constant one can construct solutions where the matter part occupies more 
than half the three sphere. Understanding these solutions requires a proper 
coordinate system covering the entire space. 

Let us now insert the constant of integration expressed as central values into 
the interior metric 

ds 2 = -(attR 2 (p c +Pc)(1 - cosx) + cosx) dt 2 



R 2 d X 2 + R 2 sm 2 xdtt ■ (3.12) 



While we have been able to find the metric, in order to find all unknown func- 
tions, we have to impose one additional condition to close the system. As this 
final condition (Condition 3) we assume q — er 2 — eR 2 sin 2 x, where e has di- 
mension charge/length 2 . This condition immediately implies that the effective 
energy density is constant, by virtue of (I2.17[) . Note that the resulting energy 
momentum tensor is regular everywhere, however, the proper charge density a 
blows up like 1/r near the center by virtue of (|2.2p . 

This allows us to find the radial pressure of the matter configuration as a 
function of the radius or the third Euler angle. For instance, using p = z — p 
we obtain 

= (Pc + Pc)(^Pc - A - e 2 ) - p c (4ir(p c + 3p c ) - A - e 2 ) cosx ,g ^ 
P (47r(/9 c + 3p c ) - A - e 2 ) cos x - 12tt(p c +p c ) 

This pressure function vanishes at x = Xb where 

(p c +p c )(^p c - k-e 2 ) 

cosx& = 7-: , -r-r 7 57 ■ (3.14) 

p c (47r(p c + 3p c ) — A — e 2 ) 

Note that x& does not have to exist. In fact, we will construct solutions where 
the pressure never vanishes. 

Since we have found an expression of the form \b = XbiPcPc, A, e 2 ), we can 
also express the central pressure in terms of the other quantities. For this we 
find 

= PcQ- - cosxb)(47rp c - A - e 2 ) 
Pc 4TT Pc {3cosxb-l)+A + e 2 ' l ' ' 

Since we are only interested in matter distributions with regular center, we 
assume the central pressure to be finite which is equivalent to the positivity of 
the denominator of (|3.15p . Therefore, we find 

^ 1 A e 2 ^ i6 ^ 

3 Ylirp c Ylitpc ' 
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which is a generalization of the well-known Buchdahl inequality. To see this, we 
go back to the radial coordinate r and use that cos(arcsinx) = yl — x 2 . This 
gives 



cos (arcsin ^) = Jl - ^ = e~ b ^' 2 , (3.17) 

by virtue of (|2 . 18|) . Our Condition 3 states that q/r 2 — e, inserting this 
into (f2T20|) leads to 



-b(r b ) _ , _ 8vrp 2 _ e 2 _ A 2 



1 - -irK - -r 2 - -r 2 . (3.18) 



3 3 3 
Let us finally denote M = 47r/3 pr^, then equation p. 161) becomes 

A e 2 
3 127rp c 127rp c ' 




(3.19) 



which is the familiar form of the Buchdahl inequality, see also [HE]. Letting 
A, e 2 — > gives the textbook result 



2M 1 2M 8 , 

1 ~>8- ir<5- (3 - 20) 

3.3 Classes of new solutions 

Based on the pressure function (|3.13l) we can straightforwardly identify new 
classes of solutions. Firstly, let us assume that the pressure vanishes at the 
equator of the three sphere, this means when Xb — tt/2. This type of solution 
occurs when the relation A + e 2 = 47rp c holds. When A + e 2 > 47rp c , the 
pressure vanishes after the equator of the three sphere. If we furthermore want 
the pressure to vanish before the second center of the three sphere, we must 
impose 47rp c < A + e 2 < 8Trp c (p c + 2p c )/(2p c + p c ). 

Next, we can encounter situation where the pressure does not vanish any- 
where and the solution has a second regular center, this means we assume the 
pressure at the second center to be finite. In that case 8irp c (p c +2p c ) / (2p c +p c ) < 
A + e 2 < 8ir(2p c + 3p c ). If A exceeds this upper limit the pressure diverges. In 
this class of solutions, there is one special solution where the pressure is constant. 
This happens when A + e 2 = 47r(p c + 3p c ) which is a charged generalization of 
the Einstein static universe. In the following Section, we will discuss these new 
types of solutions in detail. 



4 New exact solutions 

4.1 Solutions with exterior Nariai metric 

Let us assume A + e 2 = 4irp c . In this case the pressure function p.!3[) vanishes 
when \b = 7r/2, the equator of the three sphere. Now we need to identify the 
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correct vacuum solution matched at the vanishing pressure surface. It is clear 
that the vacuum solution in this case cannot be part of the Reissner-Nordstrom 
de Sitter spacetime. The only other static and spherically symmetric spacetime 
with charge and cosmological constant is the charged Nariai metric, we collect 
some basic facts of this solution in the Appendix, [S] 
This metric is given by 



i-(-(asin?A + ^cosV>) 2 di 2 +^ 2 ) + ^dQ 2 , (4.1) 



ds 

A 

where the cosmological constant and the total charge are related to A and B by 

A = \{A + B), Q 2 = l -{B-A), (4.2) 

A = A-Q 2 , B = A + Q 2 . (4.3) 

This is the most general form of the charged Nariai solution, the constants a 
and /3 can be chosen arbitrarily and will be fixed by our matching conditions. 
In the following we will show that metric (|4.ip is the correct exterior for (|3.11[) 
or (|3.12[) when the matching is performed at the \ — 7r /2 hypersurface. 

There are two equivalent approaches to matching two metrics at a given 
hypersurface. One can either match the first and second fundamental forms or 
one can introduce Gauss coordinates relative to the hypersurface and then show 
that metric is continuous and differentiable at the matching surface. We will 
follow the latter approach since both metrics are almost in the correct form. 

Let us introduce a new coordinate £ = R\ for (|3.11[) . which we will call the 
interior metric 

d*L = - ( 1 + ^ C + °y /E) ) "dt 2 + d( 2 + R 2 sin 2 ((/R)dil 2 . (4.4) 

Let us also introduce a new coordinate ( = i/j/y/A for metric (|4.1[) . the exterior 
metric 

d 4rt = -\{asm(VAQ + (3 cos{VA (jf dt 2 + d( 2 + \-dtt 2 . (4.5) 
A B 

The matching surface (vanishing pressure surface) is located at Q — wR/2. 
Continuity of the metric yields 

R 2 = I (4.6) 

= -^={asm(^TrR/2) + f3cos(y/AnR/2)) (4.7) 

while continuity of the first derivative implies 
2 C 2 



. (asm(VATrR/2) + /3cos(V^7ri?/2)) 
1 + G H A 

x (av / lcos(v / A7ri?/2) - pVAsm(VAirR/2)) . {4, 
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It is straightforward to solve these equations for a and fi. The most elegant 
approach seems to be to firstly use that R = l/\/B, and secondly use the 
second equation for the first factor of the right-hand side of the third. This 
gives 



/\ZAtt\ /VAn\„ 2C 

C °<7^2) a - S <7B ;2> = -R- (4 ' 10) 



In the absence of charge, A = B = A and these equations simplify considerably 
since the trigonometric functions become either 1 or 0. Let us denote £ = 
^A/Bn/2, then one can write the latter equations as the following simple 
linear system of equations 



sin£ cos£ \ f a\ A/I/(l + C) 
cos£ -smi) \fi) \ -2C/R 



(4.11) 



The matrix on the left is its own inverse and thus a and fi are given by 

(4.12) 



sin<e cos£ \ /VZ/(1 + C) 
fi) _ -sin£/ ^ -2C/R 

Therefore, we have shown that the interior and the exterior metric can be 
matched, the metric being C 1 at the surface. Without further assumptions, 
this cannot be improved, as can be seen by the following argument. The energy 
density inside the charged star is constant, while it is zero in the vacuum region. 
Therefore, the energy-momentum tensor has a jump. By virtue of the Einstein- 
Maxwell field equations, the Einstein tensor must have a jump too. Since it 
contains the second derivatives of the metric, the metric is at most C 1 . 

We would also like to note that these interior solutions can be interpreted 
as matter solutions which require the Bertotti-Robinson spacetime [14, .18] to be 
the electro- vacuum part of the manifold if we assume a vanishing cosmological 
constant. 



4.2 Solutions with black hole event horizons 

If 47rp c < A + e 2 < 8irp c (p c + 2p c )/(2p c + p c ), then the pressure vanishes 
after the equator of the three sphere but before its second center. In this case 
the vanishing pressure surface is in a region where the area group orbits are 
decreasing. The vacuum part of this spacetime is part of the Reissner-Nordstrom 
de Sitter solution. However, this region contains the black hole event horizon. 
This unusual class of solutions has been studied in the past in the absence of 
charge, see [5][6]. Since these solutions appear to have little physical relevance, 
we will not discuss them further. 
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4.3 Generalized Einstein static universes 

We now assume that A + e 2 is in the range 8np c (p c + 2p c )/(2p c + p c ) < A + e 2 < 
8n(2p c + 3p c ). In this case the pressure is strictly positive. The spacetime has 
two regular centers and it is natural to refer to such spacetimes as generalizations 
of the Einstein static universes. 

Within these solutions, there is the charged generalization of the original 
Einstein static universe. If we choose A + e 2 = Air(p c + 3p c ), the prefactors 
of the trigonometric functions in the pressure (|3.13l) vanish and the pressure 
is constant. The radius of this charged Einstein static universe is 1/i? 2 = 
4ir(p c +p c ). 

It should be noted that in all of the above solutions, the cosmological con- 
stant and the charge parameter appear in the combination A + e 2 and one can 
regard this as an effective cosmological constant, note that e 2 > 0. For the 
charged Einstein static universe in particular this means that we can construct 
such a matter configuration without cosmological constant. In this case the 
charge acts as an external force that pulls the matter apart and balances the 
gravitational attraction. 

Solutions with two centers are symmetric with respect to x — tt/2. To see 
this, let us denote the pressures at the two centers by pi and P 2, respectively, 
then 

Pi :=P(X = 0) =p c , (4.13) 
The latter equation can be solved for A + e 2 , which then becomes a function of 

Pc-.Pl-.P2 

2 8irp c (p c + 2 Pl + 2p 2 ) +3p!p 2 . . 

A + e = ■ ■ . (4.15) 

2/5 c + Pl + P2 

This allows us to eliminate A + e 2 from the pressure function (I3.13|) and we 
arrive at 

/ u i Pcipi +P2) +2p!p 2 + p c (pi -_p 2 )cosx 1C . 

p(x)lPc,Pi,P2\ = — : : ■ (4.16) 

2p c +P1+P2- (Pi - P2 ) cos x 

The pressure vanishes at Xb where 

cosx & = ^-P;)^- 1 )-^) (4.17) 
Pi (Po +Pi)+P2 (2p Pi -P0+P1) 

The aforementioned symmetry property of the pressure function can now 
easily be expressed mathematically, namely we have 

p(ir/2-j)[p c , Pl ,p 2 } =p(%/2 + j)[p c ,p 2 , Pl ]. (4.18) 
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5 Conclusions 



The aim of this paper was to construct new classes of exact solutions of Einstein's 
held equations in a static and spherically symmetric setting. We considered a 
charged perfect fluid to be the matter source. One class of new solutions is an in- 
terior solution which requires the charged Nariai metric to describe the vacuum 
part of this spacetime. The other class generalizes the Einstein static universe. 
These solutions are characterized by two regular centers and a non-uniform pres- 
sure function (|4.16l) . They possess an additional symmetry property namely, 
the pressure is point-symmetric with respect to the middle between the two 
centers. 

In our particular setting the cosmological constant A and the charge e always 
appear in the combination A + e 2 which we can either view as an effective charge 
or an effective cosmological constant. It is the latter point of view which might 
be of interest in the context of cosmology. If these results are not an artefact of 
our setting, then a small overall charge density throughout the universe might 
be viewed as a natural explanation of a small cosmological constant. 
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A The charged Nariai solution 

The charged Nariai solution forms a two dimensional subspace of the Reissner- 
Nordstrom de Sitter class of solutions. The inner, outer and cosmological hori- 
zons are all equivalent, sometimes called degenerate. Note that this implies 

I-™ + *-£-0, (A,, 
r r £ 3 
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the breakdown of the the (r, t) coordinate system which describes only part of 
the complete spacetime. We will follow the approach in [17] where it is noted 
that the Nariai spacetime can be viewed as a 4 dimensional submanifold of a 
flat 6 dimensional Lorentzian manifold with product structure. Let us consider 



ds = -dxg + dx\ + dx 2 + dx z + dx\ + dx% , (A.2) 



such that 



• T o + X l + X 2 — > X 3 + x \ + x \ — ~g ' (-A--3) 

with A 7^ B = const, and A + B = 2A. Thus this spacetime is a direct product 
of two 3 dimensional manifolds, one being a 3 sphere and the other one being a 
hyperbolic 3 space. To make this relation explicit, let us consider the following 
parametrization 

x Q = —= (a sin tp + (3 cos x) sinh( V a 2 + /3 2 t) , 

VA y/ a 2 + P 2 

xi = ^=^=L==(asin-0 + /3cosx) cosluVa 2 + P 2 t) , 
VA ^Jo? + p 2 



X2 = 7t ^ (/ 3sin V' - acosx) , 

VA Va 2 + (3 2 



X3 = —= sin 9 cos <j> , 
vB 

Xi = —= sin 9 sin <fi , 
\B 

x 5 = —= cos 6 . (A.4) 

Vb 



Metric (|A.2|) then gives metric (|4.1|) . 
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